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The UNSW Extrasolar Planet Search: Methods and First 
Results from a Field Centred on NGC 6633 
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ABSTRACT 

We report on the current status of the University of New South Wales Extrasolar 
Planet Search project, giving details of the methods we use to obtain millimagnitude 
precision photometry using the 0.5 m Automated Patrol Telescope. We use a novel 
observing technique to optimally broaden the PSF and thus largely eliminate photo- 
metric noise due to intra-pixel sensitivity variations on the CCD. We have observed 
8 crowded Galactic fields using this technique during 2003 and 2004. Our analysis of 
the first of these fields (centred on the open cluster NGC 6633) has yielded 49 vari- 
able stars and 4 shallow transit candidates. Follow-up observations of these candidates 
have identified them as eclipsing binary systems. We use a detailed simulation of our 
observations to estimate our sensitivity to short-period planets, and to select a new 
observing strategy to maximise the number of planets detected. 

Key words: planetary systems - binaries: eclipsing - methods: observational - meth- 
ods: data analysis - open clusters and associations: individual: NGC 6633. 



1 INTRODUCTION 

The detection of extrasolar giant planets via high-precision 
radial velocity measurements of dwarf stars in the Solar 
neighbo urhood has been highly successful. Since the first de- 
tection jMwor^^^ucloz 199 51) over fOO planets have been 
found (e.g. lButler et aJ2003tlTinnev et all2003l:lPepe et alJ 
120041) . About 20% of these planets have periods of less than 
10 days. 1 These "hot Jupiters" have a significant probability 
(~ 10%) of having an orbital inclination that causes them to 
transit their host star. This makes possible an efficient pho- 
tometric method of detecting such planets. In recent years 
a number of groups have begun photometric monitoring of 
large samples of stars, with t he aim of detecting transiting 
hot Jupiters (see Home 2003 for a list of projects). The first 
known transiting planet, HD 2094 58b, was found by the ra- 
dial velocity method |charbonneau et al.l2000l : lHenrv et alJ 
120001: iMazeh et alJl200CTl . Only in the past year have tran- 
sit searches finally begun to bear fruit, yielding five plan- 
ets from the OGLE survey iudalski et alJl2002allblcl pEl 
Konacki et all l2003l : iBouchv et alJ 120041: iPont et alJ hoo4 
Konacki et ail [2005). and th e first confirmed d etection by 
a small, wide-field telescope jAlonso et aflEooil) . 



* E-mail: mgh@phys.unsw.edu.au 

1 http://www.obspm.fr/encycl/catalog.html 



Due to geometric and time-sampling considerations, the 
transit method has a lower probability of detecting a given 
planet around a star than the radial velocity (RV) technique. 
However, a search for transiting planets has considerable 
advantages. 

• A transit search can be performed with a large-format 
CCD on a wide-field telescope, and can therefore monitor 
a large sample (~ 10 4 ) of stars simultaneously. For planets 
with periods shorter than ~ 10 days, this outweighs the low 
transit probability and makes a transit search potentially 
more efficient than a radial velocity search, which necessar- 
ily targets fewer stars. Transit searches therefore have the 
capacity to build a larger sample of known close-in giant 
planets. 

• Measurements of RV variations in the host star yield 
most of the planet's orbital parameters, but only a mini- 
mum mass (MpSini). For a transiting planet, the orbital 
inclination (i) can be estimated from the lightcurve (typi- 
cally to within a few degrees), therefore the actual mass of 
the planet can be measured. 

• The radius of the planet can be measured from a fit 
to the transit lightcurve (and an estimate of the host star's 
radius), therefore the average density of the planet is known, 
providing an important constraint for models of its structure 
and composition. 
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• Transit searches can be tailored to search for plan- 
ets around s tars in various environments, such as in ope n 
clusters (e.g. iMocheiska ct al. 2004; von Braun et al"ll2005lh 
globu lar clusters (e.g. IGilliland et al.ll2000l : IWeldrake et alJ 
2004), or in the Galactic disk. Radial velocity searches need 
to target relatively bright stars, and are therefore better 
suited to searching for planets in the Solar neighbourhood. 

• With ultra- high-precision photometry (< 0.1 mmag) , 
transiting Earth- like planets may be found ('e.g. lDeegl 2002'). 
These will be key t argets for planned space m issions such as 
ESA's Eddington jRoxburgh fc Favatall2003T) and NASA's 
Kev ler llBorucki et alJ^OO.tl. 

• iRavmond. Quinn fc Luninel l)2004h have shown that ter- 
restrial planets can form in the habitable zone (where liquid 
water can exist on the surface of a planet) of a star host- 
ing a hot Jupiter. Assuming that the orbits in a planetary 
system are near co-planar, stars with known transiting hot 
Jupiters offer an excellent opportunity to detect (or rule out 
the existence of) terrestrial planets in the habitable zone. 

• While models suggest that planets may form and 
remain in stable orbits around binary stars (e.g. 
lArtvmowicz fc Lubowl Il994l : iHolman fc WiegerU Il99flh . no 
circum-binary planets are known to date. 2 Such planets 
could be found by targeting known eclipsing binary systems, 
which offer an increased probability of a transiting configu- 
ration (if the orbits are co-planar), as well as the possibility 
of detecting n on-transiting planet s via precise timing of bi- 
nary eclipses (|Dovle fc DeedEooiL and references therein). 

• If the host star is sufficiently bright, the atmosphere of 
the planet can be studi ed using high-precision spectroscopy 
during a transit fe.g. |fji ch neider| |l994|: eager fc Sasselovl 
l2000l:lBrowrJ200ll:IWebb fc Wormleatonl200llh To date, this 
has been done successfully for the planet orbiting HP 20945 8 
ijCharbonneau et aljEoOa IVidal-Madiar et aLlEiol l2004lh 

• Infrared emission from a transiting planet may be de- 
tectable by the disappearance of s pectral features when th e 
planet is behind the host star (e.g. iRichardson et al.ll2003F) . 

• Ultra-high-precision photometry and timing of plane- 
tary transits has the pot ential to detect moons and rings 
around giant planets (e.g.lBrown et al.ll2001i iDovle fc Deed 
l2004l:lBarnes fc Fortnevll2004h . 

• At a photometric precision of 10 _J , a transiting gas 
giant's oblateness can have a detectable effect on the 
lightc urve, placing c onstraints on the planet's rotation rate 
iSeaeer fc Huil2Q02l) . Lensing of the host star 's light by the 
plane t's atmosphere may also be detectable llHui fc Seagerl 
2002), providing further information about its structure. 

Transit searches ideally require a wide-field telescope, 
in order to monitor a large sample of relatively bright stars. 
Transiting planets around stars bright enough for high- 
resolution spectroscopic follow-up are the most valuable de- 
tections. 

The other key ingredient is the ability to obtain high- 
precision photometry. A typical hot Jupiter (radius 1 Rj up , 
orbiting a Sun- like star with a 3-day period), will cause a 
1% decrease in the host star's apparent brightness for ~ 
3 hours. Therefore a transit search needs to achieve a relative 



2 With the possible exception of the system HD 202206, where 
the outer planet may have formed in a circu m-binary disk aroun d 
the primary and a brown dwarf companion iCorreia et al 12004) . 



photometric precision of a few millimagnitudes over time- 
scales of at least a few hours, which is a challenge for CCD 
photometry with a wide-field, ground-based telescope. The 
present paper is primarily concerned with the methods we 
have used to reach this precision, and the first results of our 
search. 

In section |5J we describe the Automated Patrol Tele- 
scope and the camera. In section [3] we identify the major 
factors limiting the photometric precision. The techniques 
we have used to overcome some of these limitations are 
discussed in sections 2] an d El In section [5] we present the 
data we have obtained to date, the quality of the resulting 
lightcurves, and the first results from our search. Estimates 
of our expected detection rate in the current data, as well as 
for our improved observing strategy, are given in section Q 



2 THE AUTOMATED PATROL TELESCOPE 

The Automated Patrol Telescope (APT) at Siding Spring 
Observatory, Australia, is owned and operated by the Uni- 
versity of New South Wales (UNSW). It is a 0.5 m telescope 
of Baker-Nunn design, with a 3-element correcting lens and 
an f/1 spherical primary mirror. It was originally used as 
a photographic satellite-tracking camera, before being do- 
nated to UNSW by the Smithsoni an Institution, and con- 
verted for use with a CCD camera JCarter et alJll992lh 

The camera was built by Wright Instruments, Enfield, 
U.K.. It uses an EEV CCD05-20 chip with 770 x 1150 
(22.5 [im) pixels to image a 2 x 3 degree field. However, 
the telescope has a potential field of view at least 5 degrees 
in diameter, and a new camera is being built in 2005 to take 
advantage of this. 

The telescope is housed in a building with a roll-off roof. 
The operation of the roof, the telescope and the CCD cam- 
era are all controlled via a PC running Linux. An Internet 
connection to the PC gives a remote observer full control of 
the telescope. A local weather station and several web cam- 
eras allow the observer to monitor weather conditions and 
the operation of the telescope. The system also has the abil- 
ity to close automatically in the event of rain, high humidity, 
or loss of connection to the Internet. 



3 FACTORS LIMITING PHOTOMETRIC 
PRECISION 

We have used simulated APT images to study the factors 
contributing to the photometric noise. Here we summarise 
only the most important factors for bright stars. 

3.1 Poisson noise 

For unweighted aperture photometry, the variance in flux 
due to Poisson noise is up j = (f + s)/g, where / and s are 
the counts in the effective photometry aperture due to the 
star and sky respectively (<7p,f, f, and s all in ADU), and g 
is the CCD gain (g « 8 e^/ADU for the APT). Expressed 
as an RMS variation in magnitude, provided apj <C /, this 
becomes 



a f + s 

" = 7 1 



(i) 
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where a = 2.5/ In 10 ~ 1.086, and the magnitude of the star 
is m — z — 2.5 log /. The magnitude zero point, z, is the 
magnitude of a star which results in one ADU of detected 
flux at zero airmass. For a 2.5-minute APT image in V band, 
z ~ 22.5. The precision we reach in practice is compared to 
this limit in Fig. 2] 

3.2 Scintillation 

Scintillation sets the theoretical minimum noise level for the 
brightest (unsaturated) s tars. The magnitude scatter due to 
scintillation is given by JKieldsen fc Frandsenlll992t equa- 
tion 3) 

a scint = (0.09mag)D- 2/3 x 3/2 At- 1/2 e-' l/8 (2) 

where D is the telescope diameter in centimetres, \ m the 
airmass, At is the exposure time in seconds, h is the altitude 
of the observatory in km. Using typical values for our ob- 
serving program with the APT (D = 50, x — 1-5, At = 150, 
h — 1.1), the estimated scintillation limit is 0.9 mmag RMS 
per image. 

3.3 Undersampling and intra-pixel sensitivity 
variations 

The CCD currently used in the APT camera has 9.43 arc- 
sec pixels. The image of a star in the focal plane has a 
FWHM equivalent to only ~ 0.7 pixel (Toyozumi & Ash- 
ley, in preparation). It is far from a Gaussian profile, hav- 
ing a narrow, bright core, and broad, asymmetric wings 
(Fig. [5] left panel). This "instrumental point spread func- 
tion" (iPSF) is largely due to the telescope optics, and is 
therefore stable and largely independent of local seeing con- 
ditions. Its shape does, however, vary slightly with position 
on the CCD. When sampled by the CCD pixels, it becomes 
the "effective PSF" (ePSF), with a width of ~ 1.0 pixel. For 
a focused star centred on a pixel, about 52% of the light falls 
within that pixel. 

Because of this severe undersampling, the effect of CCD 
intra -pixel sensiti vity variations (IPSV) becomes important 
(e.g. lLauerlll999l) . Direct measurements by scanning the 
APT CCD with a small (< 0.2 pixel) spot of light show that 
the point se nsitivity varies by more th an 30% over a pixel in 
the V band jTovozumi fc Ashlevl2o'o5T) . As a result, the total 
flux detected for a star varies with its precise position within 
a pixel (Fig. [31 left panel) . Combined with inevitable small 
shifts due to telescope jitter and differential refraction, this 
constitutes a significant source of photometric noise. Since 
IPSV is largely caused by the surface structure of the CCD, 
its amplitude decreases with increasing wavelength. For fo- 
cused images, the detected flux varies by 4% in V band, 
and 2.5% in / band. Due to the uniform structure of the 
CCD, the sensitivity variation is, to a good approximation, 
identical in all pixels. 

Conventionally, PSF fitting is the technique of choice for 
general purpose precision photometry. However, with under- 
sampling, IPSV, and variation of the PSF over the field of 
view, photon-noise limited PSF fitting is extremely difficult 
to achieve. In principle these effects can be accounted for, 
and practica l me thods for doing so have been described by 
iLaueJ lll999l) and lAnderson fc Kind fcOOd) . We have applied 



the algorithm of Anderson & King to APT images. While the 
the result was a significant improvement over naive aperture 
photometry , a precision below ~ 1% could not be reached for 
most stars llHidas et alJl2003l . l2004h . The severe undersam- 
pling of APT images is the likely cause of this limit. Most of 
the objects to be measured are blends of several unresolved 
stars. Thus the true number of PSFs contributing to the 
observed pixel flux distriubtion is often significantly higher 
than what the model is fitting (or is capable of deriving from 
the image itself). 

3.4 Flatfielding 

We have made sure that flatfielding errors do not make any 
significant contribution to the photometric noise for bright 
stars. Any changes in the pixel-to-pixel sensitivity variations 
over time-scales of a few months are at a level below 0.1%. 
We combine a sufficient number of twilight flatfield images 
so that the Poisson noise in the combined master flat is also 
well below 0.1% RMS per pixel. 

Note that as twilight skies have generally bluer optical 
colours than dark skies, while most target stars generally 
have redder colours, the pixel response as a function of wave- 
length also needs to be taken into account. However, this is 
not a problem if only the overall sensitivity is dependent 
on wavelength, and the relative pixel-to-pixel differences are 
constant. This seems to be the case for most CCDs. 

3.5 Non-linearity of CCD response 

The relationship between the flux incident on the CCD and 
the resulting pixel values can be significantly non-linear, 
largely due to the CCD readout amplifier. This can be an- 
other limiting factor for precision differential photometry 
from undersampled images, where the distribution of pixel 
values in the photometry aperture can vary considerably 
from image to image. 

Fig. shows our measurement of the response of the 
APT CCD to a range of incident flux levels (in flatfield im- 
ages). The signal in bright pixels deviates from linearity by 
up to 6%, which is large compared to many CCD cameras. 
We have fitted two third-order polynomials (above and be- 
low 13000 ADU) to these data, and use this fit to correct 
every image before further processing. As the residuals from 
the fit are less than 0.05%, this correction reduces the effect 
of non-linearity on the photometric noise to below 0.5 mmag. 

3.6 Crowding 

In a crowded field, photometry apertures placed around 
stars frequently have relatively bright pixels near their 
boundaries. These are from blended or nearby stars within a 
distance of ~ 30" (the radius of the aperture we use) . Small 
shifts in position (~ 1 pixel), due to telescope jitter and im- 
age rotation, cause variations in both the brightness of such 
pixels, and the extent of their overlap with the aperture. 
The amplitude of this additional noise will depend on each 
individual star's immediate neighbourhood in the field. In 
the worst case, the presence of a relatively bright neighbour 
can add 10-20 mmag (RMS) of noise. This can occur for all 
stars, but will occur more commonly for the fainter stars. 
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Figure 1. Deviation from linearity of the CCD response, mea- 
sured from a set of dome flat images taken at a range of exposure 
times. The fitted curve (two third-order polynomials joined at 
~ 13000 ADU) is used to correct images for this effect. 



Furthermore, the presence of faint, unresolved back- 
ground stars makes the effective sky level more noisy and 
variable over short spatial scales. Thus, while crowding does 
not place a strict limit on the precision, on average it in- 
creases the level of noise in the lightcurves. 

In order to minimise the additional photometric noise 
due to crowding, we have experimented with allowing each 
object to have a different sized photometry aperture. The 
aperture radius is chosen to minimise the flux in the bright- 
est pixel near the boundary in a sample image, while ap- 
plying a weight towards smaller apertures. Once a radius is 
chosen for an object, it is kept fixed for all measurements. 
For most stars fainter than V = 14.5 this reduces the RMS 
scatter by at least 20%. To date we have found no clear im- 
provement for the brighter stars, though future refinement 
of this idea may lead to significant gains. 

Besides contributing to the photometric noise, crowd- 
ing has another important effect on our ability to detect 
variations in a star's magnitude (including transits). The 
additional light from blended stars "dilutes" the amplitude 
of any variations, reducing the signal to noise ratio of plane- 
tary transits, and making the eclipses of binary stars appear 
more like transits (see section IS3I . 



3.7 Atmospheric conditions 

Relative photometry is generally not affected by temporal 
variations in overall sky transparency over the course of a 
night. However, the transparency can also vary spatially, 
for example due to undetected thin clouds. For a wide-field 
instrument like the APT, this can result in significant varia- 
tions in magnitude differences across the field, and thus ad- 
ditional noise in the relative photometry. When such "non- 
photometric" conditions occur, they can degrade the preci- 
sion for bright stars (9 < V < 13) to 5-10 mmag, making 
the data obtained only marginally useful. 



3.8 Systematic errors 

Systematic variations in apparent magnitude common to 
all stars are removed by calibrating the photometry us- 
ing a large sample of stars (see section 15.211 . However, 
in some of the resulting lightcurves, systematic trends 
do remain. These are usually in the form of a linear 
change in apparent magnitude with time, amounting to 
a ~ 10 mmag decrease or increase over the course of a 
night. We note that similar trends have been observed in 
lightcurves from oth er photomet ric surveys, such as HA TNet 
jBakos et alJ 120041,. MACH O forake fc Cookl 12004 . and 
OGLE JUdalski et alJl2003l) . though the latter authors also 
cl aim that these can be removed u sing a method developed 
bv lKruszewski fc Semeniukl i2003T) . 

In our own data, these systematics may be largely due 
to the inability of our sky fitting routine to accurately fol- 
low changes in the effective background signal over small 
spatial scales. This can occur in regions of the field where 
the density of faint, unresolved background stars is large 
(e.g. near the centre of an open cluster). Combined with the 
time-dependent slope in the intrinsic sky brightness, this can 
lead to the kind of trends we observe. We are working on an 
improved sky fitting algorithm in order to minimise these 
effects. 

Small changes in the camera focus also lead to system- 
atic errors in the relative photometry, as they alter the con- 
tribution of light to each star's photometry aperture from 
blended or nearby stars. Since the magnitude of this change 
varies depending on each star's immediate neighbourhood, 
it cannot be simply corrected for. While the focus is sta- 
ble over the course of a typical observing run, it does need 
to be reset occasionally, in particular after a filter change. 
This effect can cause magnitude offsets from night to night 
in the lightcurves of some stars (those whose fractional con- 
tribution from neighbours is significantly different from the 
reference stars used for calibration). 



4 OBSERVING STRATEGY 
4.1 The raster-scan technique 

Tests with severely defocused images (FWHM « 10 pixels) 
have shown that a relative photometric precision of at least 
~ 3 mmag was possible with the APT, for bright (V ~ 7) 
stars in sparse fields. Of course, such images are not ap- 
propriate for observing a large number of fainter stars in 
crowded fields. 

A new observing technique has allowed us to achieve 
the same precision without degrading image quality (Ashley 
et al., in preparation). During each exposure, the telescope 
is moved in a precise raster-like pattern covering the area 
of a single CCD pixel. This broadens the effective PSF only 
slightly (FWHM ss 1.5 pixels, see right panel of Fig. 
However, because each point on the instrumental PSF scans 
over a range of sub-pixel positions during the exposure, the 
sensitivity variations within the pixel are effectively aver- 
aged out (Fig. Residual variations in the detected flux 
due to IPSV are limited to < 1 mmag. 

A slightly different technique, based on the same idea 
of broadening the PSF by moving the telescope during 
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Figure 2. The instrumental PSF (iPSF) of the APT near the 
centre of the field. Each PSF has been normalised to have a total 
flux of 1. The contour levels are identical in both panels, and have 
uniform logarithmic spacing. The grid spacing corresponds to 1 
CCD pixel. Left panel: The normal, focused iPSF, determined 
by fitting a multi-Gaussian model to a set of 150-sec, V-band 
images of a bright star (Toyozumi & Ashley, in preparation). The 
FWHM is ~ 0.7 pixel and the peak intensity ~ 1.1. Right panel: 
The broadened iPSF obtained by the raster-scan technique. This 
model was calculated by numerically convolving the normal iPSF 
with the raster-scan pattern followed by the telescope (Ashley et 
al., in preparation). It has a FWHM of ~ 1.5 pixels, and a peak 
intensity ~ 0.5. 




-0.04 -0.02 0.02 0.04 -0.001-0.0005 0.00050.001 



(Magnitude) (Magnitude) 

Figure 3. The effect of intra-pixel sensitivity variations (IPSV) 
on the photometry from undersampled images taken by the APT. 
Each image represents one CCD pixel, and shows the relative 
magnitude measured for a star, as a function of its position within 
the pixel. An accurate model of the instrumental PSF, obtained 
from real images (Toyozumi & Ashley, in preparation), was con- 
volved with the pixel sensi tivity, as measured by sc anning the 
CCD with a spot of light iTovozumi fc Ashlevl [20051) . The two 
panels correspond to the two iPSF models presented in Fig. |2] 
Left panel: the result for normal, focused APT images in the V 
band. Right panel: the reduced effect when the PSF is broadened 
by the raster-scan technique ("section 14.11 . Note the difference in 
scale. 

exposure is used b y the Hungarian Automated Telescope 
l|Bakos et, alJl2004h . 

4.2 Nightly routine 

Most of the observations are conducted remotely via the In- 
ternet. The observer's main task is to initiate the observing 
process, and monitor weather conditions and the quality of 
the data obtained. Each of the tasks of starting up the tele- 
scope system, opening the roof, taking twilight flats, and the 



observations themselves are executed via simple shell scripts 
on the telescope-control PC. 

Twilight fiatfield images are taken on most clear 
evenings when the Moon is down, and on some mornings. 
~ 15 images with sufficiently high background flux can be 
taken per twilight. Other calibration images are not neces- 
sary for the APT camera. The bias level shows no signifi- 
cant structure over the CCD, and can be accurately mea- 
sured from the overscan region. Dark current is negligible 
(~ 0.6 e~ per pixel per second). 

During the night, as each image is read out, an a strom- 
etry program locates Tycho-2 stars <Hflg et aUl200ch in the 
image, and calculates the (linear) transformation from CCD 
pixel coordinates to celestial coordinates. The coefficients 
of the transformation are written to the FITS header as a 
World Coordinate System (WCS). Typically, ~ 200 Tycho- 
2 stars are matched. In this way, telescope pointing is kept 
accurate to within ~ 0.2 of a pixel. This helps reduce the 
effects of any flatfielding errors, as stars near the centre of 
the field are always imaged by the same pixels. As differen- 
tial refraction by the atmosphere causes a gradual rotation 
(~ 4 arcmin in 8 hours) and change in image scale (< 0.1%), 
stars near the corners of the field will shift by ~ 1 pixel. 

The astrometry program also performs a quick measure- 
ment of the magnitudes of the Tycho stars and compares 
them to the catalogue values. This allows the remote ob- 
server to monitor sky transparency and detect the presence 
of clouds. 

For most of our 2003-2004 observations, the telescope 
continuously cycled between four adjacent fields, taking 150- 
sec exposures, as long as the fields were observable. Includ- 
ing CCD readout and telescope slew time, each image took 
3.5 minutes. On a full clear night ~ 30 images per field were 
obtained. Due to physical limitations of the mount, the tele- 
scope can follow a field for a maximum of ~ 8 hours. 



5 ANALYSIS STRATEGY 
5.1 Data reduction 

Optimal aperture photometry can be co mpetitive with PSF 
fitting for moderately crowded fields (e.g.lNavlor 1998)). For 
brighter stars, variance-weighted PSF fi tting is equ ivalent 
to standard aperture photometry (e.g. Ilrwinl Il997l) . This 
and the constancy of the PSF for the APT have enabled 
us to develop a simple, robust and extremely accurate aper- 
ture photometry package for processing APT images. Image 
crowding, usually the bugbear of aperture photometry, is 
dealt with in a novel, yet simple manner. 

Briefly, each APT frame is processed in the standard 
manner up to and including generation of object catalogues 
using the development too lkit of the INT WFS pipeline soft- 
ware lllrwin fc Lewisj200ll) . Particular care is taken over flat- 
fielding, since we need to ensure the gross pixel-to-pixel sen- 
sitivity variations are mapped to better than the 0.1% level. 
We combine all the twilight flats (of order 100) obtained dur- 
ing a 3-4 month season. Each twilight image has ~ 20000 
counts (the gain is ~ 8 e~/ADU). 

Non-linearity of pixel counts is dealt with by first mea- 
suring the global non-linearity from a varying time series of, 
essentially, fiatfield frames under constant illumination and 
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then applying the derived constant look-up-table correction 
to all pixels before processing. 

The final stages of the INT WFS pipeline involves 
generating object catalogues for each frame. Among other 
things the object catalogues are used to derive accurate 6- 
parameter linear frame-to-frame transforms with respect to 
a chosen master reference frame. To date, a single, low- 
airmass image from a dark, photometric night has been cho- 
sen as the master frame for each field. In the future a sep- 
arate master image will be generated by stacking a dozen 
or so consecutive images. This will remove spurious objects 
from the master catalogue (such as cosmic rays), and im- 
prove the astrometric precision. 

The transformations are of the general linear form 

x = ax + by + c ; y = dy + ex + / (3) 

where a and d define the relative x and y scales, b and e 
the rotation and/or shear, and c and / the global shifts. 
Note that we could perform these transformations using 
the WCS data stored in each image header. However, this 
would not make optimal use of the information in each frame 
(the WCS transformations are calculated using ~ 200 Ty- 
cho stars, while the full object catalogues typically contain 
10000 objects to be matched to the master catalogue), and 
would limit the accuracy to ~ 0.1 pixel. 

Although we attempt to place the exposures of each tar- 
get on comparable field centres within each night and from 
night to night by performing astrometry on each image as 
it is read out, small differential adjustments of order 1 pixel 
are still required. Even with apparently perfectly aligned 
images, small temperature changes throughout nights, dif- 
ferential refraction variations as a function of airmass and 
differential field rotation all contribute significant pixel level 
distortions. However, the approximate global alignment does 
ensure that non-linear terms are negligible at this pixel scale. 

The reference frame catalogue also defines the master 
coordinate list to be used for the later analysis, and this 
is transformed (to within 0.01 of a pixel) onto each frame 
to be processed. Without this stage relative accuracies of 
~ 0.1 of a pixel are still possible, from the individual im- 
age catalogues. However, with these undersampled images, 
coordinate accuracies of 0.1 pixel are not sufficient to give 
mmag photometry for the brighter objects. This is partic- 
ularly important in crowded fields, where bright pixels on 
the boundary of the photometry aperture are frequent (see 
l3~5l. Variat ions in the contribution to the aperture sum 
from these pixels due to errors in aperture centring need 
to be minimised. Tests have shown that errors of ~ 0.1 pixel 
in the aperture positions increase the photometric noise by 
50% for the average bright (V < 14) star in a crowded field 
(compared to the near-photon-limited results we otherwise 
achieve) . 

The final ingredient is a very stable and robust 2D back- 
ground estimator, which is used to estimate and remove the 
global background in a highly repeatable manner in each 
frame. With this kit of parts each dataset can be analysed 
in (almost) exactly the same way, using the same apertures 
accurately co-located on each frame. Each aperture sum in- 
cludes sub-pixel partitioning of flux and the whole package 
runs completely automatically. 

The main advantage of this approach is that system- 
atics from crowding are the same in each frame because of 



the PSF stability and the differential coordinate precision, 
and even quite crowded fields can be successfully analysed 
since essentially identical image confusion arises in all cases. 
With this approach we are achieving differential photometric 
precision of better than 2 mmag for bright objects and are 
close to the theoretical photon noise limit at all magnitudes 

(Fig. a. 

5.2 Photometric calibration 

Global image-to-image magnitude variations are removed 
by ensemble relative photometry. An approximate magni- 
tude zero-point (good to ~ 0.1 mag) is set by comparison 
with the magnitudes of 100-200 stars from the Tycho-2 cat- 
alogue. Stars in the range 9 < V < 11 are then used for the 
calibration. There are typically ~ 1000 of these in a dense 
Galactic field. Magnitude residuals (Am) from the nightly 
mean value for each reference star are calculated, and for 
each image, fitted by a function of the form 

Am = a + bx + cy + dxy + ey 2 (4) 

where x,y are the star's pixel coordinates on the CCD, and 
a—e are constants for an image. A second order term in y 
was necessary because the CCD is larger along the y-axis, 
which corresponds to the right-ascension direction (along 
which the airmass varies more). Adding an x 2 term to the 
fit does not significantly improve the photometric precision. 
Therefore, to make the fit more robust, this term is not 
included. As this function provides a good approximation to 
the variation in airmass across the field, a separate airmass 
term is not necessary. 

This fit is then subtracted from each star's magnitude. 
Mean magnitudes and residuals and recalculated. Based on 
the RMS residuals of each star over the night's data, variable 
and badly measured stars are removed from the reference 
list, and the whole process is repeated. In most cases, the 
coefficients of the fit are effectively zero after 4-7 iterations. 
Any images for which the median RMS of stars with 9 < 
V < 12.5 remains higher than 15 mmag after calibration 
are rejected, and the calibration is repeated without them. 

The above process is applied to each night's data sep- 
arately. Any remaining night-to-night magnitude zero-point 
offsets are removed at the time the final lightcurves are com- 
piled. 

5.3 The pipeline 

The entire data reduction process, from raw images to cal- 
ibrated magnitudes, can be run automatically using a Perl 
script. The only interactive part of the procedure is the se- 
lection and combining of a suitable set of twilight flats to 
create a master Hatfield image. This needs to be done once 
for each 3-4 month season of data. Running on a Sun Enter- 
prise 4500 workstation (with 8 CPUs and 4 Gb of memory), 
the pipeline will then process a full night's data in less than 
2 hours, and several of these processes can be run in parallel. 

From each dataset, two sets of magnitudes are gener- 
ated, measured with photometry apertures 2 and 3 pixels 
in radius. The larger aperture gives a slightly lower average 
RMS for most bright stars, but is more susceptible to con- 
tamination of lightcurves by close neighbours. When a small- 
amplitude signal, such as a transit, is found in a lightcurve, 
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a comparison of the signal amplitude measured by the two 
apertures provides a quick way of identifying such contami- 
nation. 

5.4 Candidate selection 

Initially, we visually inspect the lightcurves for the brightest 
2-3000 stars in each field (down to V fa 13). This procedure 
has the advantage of being less susceptible to spurious de- 
tections due to variable stars and systematic errors than 
detection software. 

We also use the tran sit detection algorithm described 
in lAigrain fc Irwml (12004) . Lightcurves are first searched for 
individual, box-shaped transit events of a given trial dura- 
tion, and the signal-to-noise ratio (S/N) of a transit at each 
trial epoch is stored. These S/N values are then used to per- 
form a period search. This is repeated for a number of trial 
transit durations, and the combination of period, duration 
and starting epoch resulting in the highest total S/N (over 
all observed transits) is selected for each lightcurve. A num- 
ber of selection criteria are applied, based on the individual 
transit and total S/N, as well as the distribution of S/N 
values (over all trial parameters) for each lightcurve. 

To date, the effectiveness of this detection software has 
been limited by the presence of systematic trends in our 
lightcurves (se e sect ion 13.81 , which lead to f requent spuri- 
ous detections. iKovacs. Bakos fc Novesl J2004T) have recently 
published details of a trend filtering algorithm which sub- 
stantially improv es the quality of lig htcurves obtained by the 
HATNet project feakos et alj 2004'l. Our initial experiments 
applying this algorithm to APT data had promising results, 
though our implementation is yet to be refined. At present 
it does greatly reduce the amplitude of systematics while 
preserving most transit signals, but also increases the white 
noise somewhat. We have not yet fully characterised the sen- 
sitivity of our detection software, but preliminary tests have 
shown that in the absence of systematics, 10 mmag transits 
are detected in a high percentage of cases (provided multiple 
transits are present and the total S/N is sufficiently high). 

5.5 Eliminating false positives 

iBrownl (l2003t) has estimated that up to ~90% of planet can- 
didates selected by a wide-field transit search are binary 
stars. Grazing eclipses of binary stars, and eclipsing binaries 
with their light diluted by a third blended star can mimic 
the depth, and (to some extent) the shape of a planet tran- 
sit. It is essential to eliminate these false positives at the 
earliest possible stage of the follow-up process. Our follow- 
up strategy is outlined below. At each stage we focus on the 
most promising planet candidates. Depending on the avail- 
able observing time we also follow up likely binary systems. 

(i) Inspection of a Digitised Sky Survey image, to re- 
veal how much blending is present in the APT image (our 
photometry apertures are effectively ~1 arcmin in diameter, 
which is typical for wide- field transit search projects). 

3 The Digitised Sky Survey was produced at the Space Telescope 
Science Institute under U.S. Government grant NAG W-2166, 
based on photographic data obtained using the Oschin Schmidt 
and UK Schmidt telescopes. 



Table 1. Fields targeted in the UNSW planet search. 



Field name 


R.A. 


Dec. 


6 


N*(V < 14) 


-^nights 


NGC 3532 


ll h 07' 


-58°42' 


1° 


7000 


61 


Gl 


ll h 31' 


-58°42' 


3° 


7000 


65 


G2 


ll h 31' 


-56°30' 


5° 


7000 


60 


G3 


ll h 07' 


-56° 30' 


4° 


7000 


59 


NGC 6633 


18 h 28' 


+6°35' 


8° 


5000 


79 


HI 


18 h 28' 


+4°30' 


7° 


4000 


47 


HI 


18 h 41' 


+4°30' 


4° 


3000 


30 


H5 


18 h 41' 


+6°35' 


5° 


5000 


30 



(ii) Estimation of the physical parameters of the sys- 
tem, in particular the radius of the transiting object, from 
the lightcurve. The lightcurve is parametrised by the pe- 
riod (P), transit depth (AF), and the durations of the 
complete transit and the "flat bottom" part (tr and tp). 
From these, the primary star's density can be calculated 
iSeaeer fc Mallen-Ornelaal2003tl . Assuming a stellar mass- 
radius relation then yields the primary's radius (Ri), and 
the radius of the transiting object R2 = R1VA.F 1 . Assum- 
ing the transit is in front of the brightest star in the APT 
aperture and adjusting the transit depth accordingly gives 
a lower limit on the size of the transiting object. 

(iii) Photometric monitoring of the target area at higher 
spatial resolution, during a predicted transit, to identify the 
host of the transiting object, and measure the true (or at 
least less diluted) transit depth. This is done in multiple 
colours so that the colour signature of a grazing or blended 
eclipsing binary can be detected. 

(iv) Medium-resolution spectroscopic follow-up to iden- 
tify the host star's spectral class, allowing a more precise 
estimate of its radius, and therefore the radius of the planet. 
The presence of double absorption (or emission) lines, fea- 
tures from different spectral types, and in particular large 
radial velocity variations in phase with the transit signal, 
would indicate a binary or blend. 

(v) For the best candidates, high-precision radial veloc- 
ity measurements (using echelle spectroscopy) will allow the 
measurement of the transiting object's mass and confirm or 
rule out its planetary nature. 



6 RESULTS 

6.1 Fields observed 

We have targeted two groups of four adjacent fields, one 
group in Centaurus, the other in Ophiuchus. Both groups are 
within 10 degrees of the Galactic plane. The field centres are 
shown in Tabled We have centred one field in each group 
on an open cluster (NGC 3532 and NGC 6633). However, 
as each cl uster contains only a few hund red relatively bright 
stars (e.g. lKoeibloedll959t Ijeffr7eslll997ll . and many of these 
are blended or saturated, only a small fraction of the stars 
we measure are cluster members. 

We began observing most of these fields in 2003, and we 
now have two full seasons' worth of data for each of them. 
The first observations using the method described in sec- 
tion0]were taken in August 2002, targeting th e open cluster 
NGC 6633 as part of a multi-site campaign llMartrn et alJ 
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Table 3. Transit candidates in the NGC 6633 field: parameters 
obtained from follow-up spectroscopy. The error on the primary 
mass (Mi) was assumed to be 10%. The velocity amplitude of 
TR-2 was estimated from the separation of double lines in the 
single spectrum we obtained. As it is a binary system with nearly 
identical components the orbital period is twice the value mea- 
sured from the APT lightcurves. We assumed the two stars have 
equal mass. The orbit of TR-4 appears to be eccentric, but we 
have insufficient data for a rigorous Keplerian fit. Assuming a 
circular orbit gives the result quoted here. 



Star ID 


Sp Type 


Mi 


K 


M 2 






(M ) 


(km s _1 ) 


(M ) 


UNSW-TR-1 


A2V 


2.5 


27 ±3 


0.34 ± 0.05 


UNSW-TR-2 


K7V 


= M 2 


87 ±3 


0.58 ± 0.05 


UNSW-TR-3 


F7V 


1.2 


53 ± 10 


0.4 ± 0.1 


UNSW-TR-4 


F5V 


1.3 


~ 20 


~ 0.2 



2004). A total of 690 images were taken through a John- 
son V filter on eight consecutive nights (7-14 August). In 
2003 we added the three adjacent fields, and obtained an- 
other 180 images of NGC 6633 in V band (June), and 380 
in R band (July- August). During the 2004 season we con- 
tinued to observe the field in V band (March-July, 1800 im- 
ages), including some nights where we again observed this 
field only, to obtain higher- cadence measurements. About 
1/3 of the images taken in 2004 (typically those taken in 
the first part of the night) were affected by a technical prob- 
lem which has rendered them largely unusable. Finally in 
August 2004 we began observing in / band, and took 230 
images of NGC 6633. For these latter observations, the ex- 
posure time was reduced to 60 seconds, as the CCD is more 
sensitive in this band. 

Due to its Northern declination, the NGC 6633 field 
can be observed for a maximum of 7 hours at relatively low 
airmass. From experience, useful images can be taken up to 
an airmass of ~ 2.2. The minimum airmass reached by this 
field is 1.3. 

With a threshold of 4cr above the background, the 
pipeline software typically detects ~ 9000 objects in the 
field, the faintest of these having V ~ 15.5. In each image, 
the positions of ~ 7000 of these objects are matched to the 
master catalogue to obtain the coordinate transformations 
into that image, with mean residual errors of ~ 0.05 pixel. 

The relative photometric precision achieved on a single 
clear night, and over a week of observations, is shown in 
Fig. |1] For stars in the range 9 < V < 11, we routinely 
achieve 2 mmag precision. For this field, the nightly RMS is 
under 10 mmag for 2300 stars. 



6.2 Transit candidates 

The lightcurves for the brightest ~ 2000 stars (down to 
V w 13) from the NGC 6633 field were searched for 
transit-like events by visual inspection. Through this pro- 
cess we identified 4 lightcurves with transits shallow enough 
(< 100 mmag) to be considered as potential planet candi- 
dates. Their lightcurves are shown in Fig. 

Table summarises the properties of the candidates, 
including an estimate of the transiting object's minimum 
radius (see item (ii) under section 15. 5H . With relatively 
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V Magnitude 

Figure 4. Photometric precision from one night (top panel), 
and 5 nights (bottom panel) of APT observations of the field 
centred on the open cluster NGC 6633 (August 2002). The solid 
line is the theoretical estimate, including Poisson noise in the 
star flux (dot dashed) and sky flux (dashed), and a magnitude- 
independent term (dotted). The latter term (1.2 mmag) repre- 
sents scintillation, as well as any residual errors due to flatfielding, 
non-linearity, and IPSV. Stars brighter than V ~ 9 are saturated 
in a 150-second APT image. 



large radii, our candidates seemed unlikely to be planets. 
The largest extrasolar pla net known to date (HP 209458b) 
has a radius of 1.35 R JuD dBrown et al]l200ll) . However, as 
iDrake fc Cookl <|2004|) point out, there are currently too few 
transiting objects with known mass to place a firm upper 
limit on the radii of extrasolar planets and brown dwarfs. 
Though not our primary targets, brown dwarfs and low-mass 
stellar objects may also yield interesting results. We also 
note that in the V band, neglecting the effects of limb dark- 
ening (as we have done for these calculations) results in R v 
being overestimated by up to 50% llSeaeer fc Mallen-Ornelaa 
2003, Fig. 12). Therefore we deemed these candidates wor- 
thy of follow-up. 

An inspection of Digitised Sky Survey images revealed 
that each of the candidates were affected by some level of 
blending in the APT photometry aperture. 

We obtained higher spatial resolution images using the 
40- inch telescope at Siding Spring Observatory (SSO). We 
aimed to monitor each candidate continuously during a pre- 
dicted transit, alternating between the V and I bands. To 
improve the time sampling we only read out one quar- 
ter of the full CCD, corresponding to a 10 arcmin square 
field (sampled by lkxlk, (f.'6 pixels). Due to a run of bad 
weather, we only observed one of the candidates (UNSW- 
TR-2) during a transit. In this case the source of the transit 
signal is in fact the second bright star in the APT aper- 
ture (~0.8 mag fainter than the brightest one, and ~ 10" 
from it). The undiluted lightcurve (Fig.|SJ shows a V-shaped 
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Table 2. Transit candidates in the NGC 6633 field: coordinates and parameters estimated from APT data. T c is the time at the 
centre of the first transit we detected (JD-2450000). Coordinates are from the WCS solution in the APT images, and represent the 
centroid (to within ~ 1") of the object, which is generally a blend. Similarly, the given V magnitudes, with errors of ~ 0.1 mag, 
are for the blend as measured by the APT. Errors in the other measured parameters are 5 X 10~ 5 d for the period, 0.01 d for 
T c , 5 mmag for the transit depth (AF), and ~ 0.5 hr for tj- and tp (the duration of the complete transit, and of the "flat 
bottom" part, respectively). Es timates for the radius of the star (R*) and the planet (R p ) were derived using the equations of 
ISeager fc Mallen-Orneiasl 120031) . UNSW-TR-1 is a member of the cluster (NGC 6633 141). 



ID 


R.A. 


Dec. 


V 


P 


T c 


AF 


t T 


t F 


Ri 


R 2 




(J2000) 


(J2000) 




(days) 


(days) 


(mmag) 


(hr) 


(hr) 


(R ) 


(Rjup) 


UNSW-TR-1 


18 29 04.4 


+06 26 54 


11.0 


2.38079 


2496.93 


30 


3.4 


2.0 


1.4 


2.9 


UNSW-TR-2 


18 30 51.9 


+07 09 20 


12.0 


1.05837 


2462.99 


25 


1.9 


0.2 


1.5 


2.4 


UNSW-TR-3 


18 31 00.7 


+07 08 25 


13.0 


1.81749 


2499.07 


70 


2.0 


0.5 


0.8 


2.0 


UNSW-TR-4 


18 32 22.3 


+06 37 13 


12.9 


4.39995 


2465.13 


50 


3.8 


0.7 


1.2 


3.8 
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Figure 5. Phase-folded lightcurves of the 4 transit candidates 
found in the NGC 6633 field. The lightcurves include all V band 
data of reasonable photometric quality (nightly rms < 10 mmag 
and at least 10 points per night) we obtained in 2002-2004. Note 
that UNSW-TR-2 is folded at its originally estimated period. 
Follow-up observations have shown that this is a binary system 
with two nearly identical eclipses, and the true period is twice 
this value. 



eclipse with a depth of 0.14 mag, suggesting that the signal 
is from a grazing eclipsing binary system. 

Simultaneously with the high-resolution imaging, we 
obtained spectra using the Double-Beam Spectrograph 
(DBS) on the 2.3 m telescope, also at SSO. We used the 




3171 3171.05 3171.1 

JD - 24 50000 



Figure 6. The undiluted eclipse signal from transit candidate 
UNSW-TR-2, observed at high spatial resolution with the 40- inch 
telescope at SSO. We observed a single transit, and alternated 
between V band (triangles) and / band (squares). The solid line 
is the model lightcurve for an eclipsing system consisting of 2 
K7V stars (with masses 0.58 Mq and radii 0.64 M Q ) in a 2.12- 
day orbit, at an inclination of ~ 83°. The dashed lines show the 
result of varying the radius of the stars by ±5% (keeping other 
parameters constant). 



highest resolution grating in both arms. The approximate 
spectral ranges covered were 5700-6700 A in the red arm 
and 3900-4400 A in the blue, with resolutions of 1.1 A 
(45 km s" 1 ) and 0.5 A (40 km s _1 ) respectively. We mea- 
sured the radial velocities by cross-correlating the H a ab- 
sorption line in each spectrum with a spectrum of the RV 
standard star (5 Vir (obtained with the same instrument). 

We also obtained a single low-resolution (2.2 A) spec- 
trum of each object covering the range 3420-5350 A. These 
were visuall y compared to templates from the UVILIB spec- 
tral library (Pickles 1998), yielding a spectral type estimate 
for each candidate (Table 

Due to the limited observing time, we were unable to 
monitor all stars blended within each APT photometry aper- 
ture. Therefore we observed the brightest star in each aper- 
ture at several independent phases, and took single spectra 
of the fainter neighbouring stars. If the transits were in front 
of the fainter stars, they would be more likely to be of stellar 
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origin (since the undiluted transit depth would have to be 
larger) . 

For three of the candidates (UNSW-TR-1, 3 and 4), the 
brightest star shows large RV variations in phase with the 
transit signal. The velocity semi-amplitude (K) of the can- 
didates is shown in Table |3] The companion masses were 
calculated using an estimate of the primary's mass from its 
spectral type, and assuming circular orbits (which is consis- 
tent with the radial velocity curve for two of these cases). 
Clearly the companions are low-mass stars, and not planets. 

6.2.1 A K7Ve binary system 

For the fourth candidate (UNSW-TR-2), the radial veloc- 
ity of the bright star dominating the light in the aperture 
showed no coherent variations when folded at the period 
measured from the transits. We only had one DBS spec- 
trum of this star's close neighbour, which was revealed in 
the 40-inch images to be the source of the transit signal. 
This spectrum is well matched by a K7V template, with 
two important exceptions: (1) there is a strong H a emission 
peak; and (2) the absorption lines and the emission line all 
have two components of approximately equal strength. Tak- 
ing into account the orbital phase at the time the spectrum 
was taken, the maximum velocity separation between the 
components is 175 ± 6 km s _1 . 

We were able to obtain two more spectra of this star at 
lower resolution. Although the double line structure is not 
completely resolved in these spectra, they do confirm that 
the separation between the lines varies with orbital phase. 
We have not performed a rigorous fit to the lightcurves. 
However, a model consisting of two identical K7V stars, each 
of mass 0.58±0.05 Mq and radius 0.64±0.01 Rq, appears to 
fit all the observations. The period is 2. 11674± 0.00002 days 
(twice the original estimate), giving an orbital radius of ~ 
0.03 AU. The orbital inclination is 83.4 ±0.1 degrees. 4 

This system may be a valuable addition to the small 
number of known double-lined eclipsing binary systems with 
late-type components, for which th e physical parameters 
can be measured with precision (e.g. iTorres fc Ribasll2002t 
iMaceroni fc Montalbanl2004l : Maceroni, C, private commu- 
nication) . Such measurements provide important constraints 
for models of low-mass stars. 

6.3 Variable stars 

Visual inspection of the individual V-band lightcurves re- 
vealed the presence of 49 well-defined variable stars. Of 
these, 3 were previously known variable stars (HD 170451, 
BP Ser and BI Oph) . We classified the remaining 46 follow- 
ing a simple procedure based on the period and the shape 
of the lightcurve. To minimise personal biases, we put stars 
into broad classes like short period pulsators or eclipsing bi- 
naries. For a few stars, we were able to determine the GCVS 
type unambiguously. 

The following types of variables were found: 40 pulsat- 
ing stars (33 short period, 5 RR Lyrae-type, 2 long period 
semi-regular) and 9 eclipsing binaries. We found a number 
of bright short period pulsating stars having A-F spectral 

4 The fitted radius and inclination are strongly correlated. 



types 5 , which suggests that they are likely 8 Scuti stars. For 
the majority, however, the lack of spectroscopic observations 
prevent more detailed classification. 

The period det ermination was performed using Pe- 
riod98 of lSperll Jl998f) and the phase dispersion minimisation 
method iStellingwerdlig?!^ . In addition to the frequencies, 
we also calculated signal-to-noise ratios in the Fourier spec- 
tra to get a me asure of reliability o f the determined frequen- 
cies. Following iBreeer et alJ ljl993ii . we adopted S/N > 4 as 
the limit of significance. The detected variable stars and 
their basic parameters (catalogue number, coordinates, pe- 
riod, type of variability) are listed in Table 0] Although in 
some cases, there is evidence for multiple periodicity, this is 
beyond the scope of the present paper. 



7 EXPECTED PLANET DETECTION RATE 
7.1 A rough estimate 

The number of detectable planets in our search can be ap- 
proximated as 

Adct = A r dwarf X Ppl X P tr X P obs (5) 

where Adwarf is the number of dwarf stars in our fields for 
which a hot Jupiter would cause at least a 10 mmag transit, 
Ppi is the fraction of stars that actually host such a planet, 
P tr is the geometric probability that the planet does transit 
(i.e. that it crosses our line of sight to the star), and P Q b s is 
the probability that we observe multiple transits, with some 
minimum combined signal-to-noise ratio ((S/N) m i n ~ 10). 

A typical transiting hot Jupiter, with radius 1.2 Rj up 
fe.g. ISozzetti et alJl20ol Fig. 3), will cause at least a 10- 
mmag transit in front of any star with radius R, 1.2 Rq, 
i.e. of spectral type ~ F8 or later. For a magnitude- limited 
sample, ~ 50% of stars are of type GKM (ICoxll2000l) . Thus 
we have Adwarf ~ 0.5 x A,, where A* is the total number of 
stars monitored with at least 10-mmag precision (i.e. with 
V < 13). For the NGC 6633 field, A* « 2000. 

The frequency of hot Jup i ters is P p i « 1% (e.g. 
iGaudi. Seager fc Mallen-Ornelasl I2005F) . and their transit 
probability is P tr ~ R*/a ~ 0.1 (a is the orbital radius). 
Therefore, observing a rich Galactic field like NGC 6633, we 
have Adet ~ Pobs- For our complete NGC 6633 lightcurves, 
Pobs ~ 0.4 (minimum one complete and one partial transit, 
1 < P < 9 days, see Fig.Q, giving an estimate of 0.4 planets 
for this dataset, or ~ 0.2 planets per year. 

In this simple estimate, we have made some gross sim- 
plifications, including 

(i) neglecting the dependence of P p i, P tr and P b s on the 
orbital period, and other properties of the star and planet; 

(ii) assuming a very simple (and perhaps optimistic) de- 
pendence of the detection probability on the transit depth, 
and neglecting the distribution of transit depths among tar- 
get systems; 

(iii) neglecting the 3-dimensional distribution of stars in 
the Galactic disk; 

(iv) neglecting interstellar extinction; 

5 Spectral types were obtained from the SIMBAD online cata- 
logue. 
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Table 4. Variable stars found in the NGC 6633 field. As this is a relatively crowded field, most objects are blends. As 
a consequence, (1) the actual variable star may be another star within ~ 30" of the catalogue star we identify; (2) the 
coordinates shown represent the centroid of the blend; (3) the V magnitudes measure the total flux in the APT photometry 
aperture; (4) the amplitudes may be diluted by blended stars, thus the values shown are only approximate lower limits. 
Coordinates were calculated from the WCS solution in the APT images, with a precision of ~ 1". The magnitude zero-point 
calibration is accurate to ~ 0.1 mag. 



Id 


Catalogue Id 




R.A. 


Dec. 


V 


P 


A 


Type 






(J2000) 


(J2000) 




(d) 


(mmag) 




VI 




18 


22 


10.2 


4-6 23 20 


12.9 


09225 


g 


G±lV_JlLr ULl 1UL1 UUlDCllUl 


V2 




18 


22 


08.3 


+6 42 42 


13.2 


78357 

KJ . 1 UOv 1 


285 


W TIMa 

V V VJ L\± Oj 


V3 


CSC 00445-00528 


18 


23 


11.7 


4-6 18 53 

L' 1 KJ ' J - J 


10.1 


0.13704 


10 


short period pulscitor 


VI 




18 


23 


06.9 


+6 42 13 


12.3 


0.14254 


10 


t;Vi(~iT'l" r"»r>T"iriH "niiltjti + rM" 


V5 


CSC 00445-00689 


18 


23 


09.1 


4-6 51 35 

I LJ ' ' 1 ' ' ' ' 


10.6 


19818 

KJ . J- rj o _L o 


9 


short period pu.ls3.tor 


V6 




18 


23 


27.5 


4-6 1 2 06 

1 LJ KJ\J 


12.8 


09452 


12 


011V_J1L [JLl 1UU \J Ll 1 o Cl L W 1 


V7 




18 


23 


33.1 


4-6 33 29 


12.8 


38673 

KJ . iJUU 1 ■ ) 


242 


W TIMa 

V V VJ Lv± cl 


V8 




18 


23 


47.9 


4-7 28 05 

I 1 — l ~ y ' ' ' 


12.4 


1.02289 


312 


Algol 


V9 


OSC 00445-01888 

v_TkJV_y IWitU U1UO O 


18 


23 


42.5 


+6 24 09 


11.1 


09855 

KJ . KJ i? O KJ <J 


10 


□ llUl l ULl 1UU l_/ Ll 1 o d L W 1 


V10 




18 


24 


27.1 


+6 45 42 


12.1 


0.19167 


13 


short period pulsfitor 


Vll 




18 


24 


32.6 


+7 30 46 


12.0 


06285 


5 


short period pulscitor 


V12 




18 


24 


40.2 


4-6 10 05 

KJ 1 V / L ' '.J 


12.5 


06559 


42 


short period pulscitor 


V13 




18 


24 


45.2 


4-6 05 31 

L' L ' (J 1 } L 


12.4 


28535 

L ' . L' '>>>'.' 


153 


RR, LyrS/C 


V14 


NOT: 6633 8 

IN VJv UUOJ o 


18 


24 


40.6 


4-7 04 05 


9.2 


1 5436 

KJ . ±Ut:iJU 


38 


□ 11U1 L UL11UU \J LI 1 o cX Lr w 1 


V15 




18 


25 


14.8 


4-6 33 54 

KJ • J- J ' ) i 


11.5 


07567 

KJ . KJ 1 t/U i 


6 


short period pulscitor 


V16 


NOC 6633 14 

1\ V_T \_J KJKJ'JiJ ± '-I 


18 


25 


04.1 


4-6 25 55 

\^ KJ iitj <JiJ 


10.0 


03658 

KJ . VJiJUiJO 


14 


clini't t~il>ti lili ni i Iqp i~rw 

JLIKJL L LJL^l 1UL1 LJ UloCLUWl 


V17 




18 


25 


10.1 


+6 46 44 


11.6 






pr"l i t^cii ti cr rii n a v\r 

K^K^Ll LJolllti L/l 11 (XI y 


V18 




18 


26 


14.4 


4-5 34 54 


13.1 


1 63336 


185 


Algol 


V19 




18 


26 


34.7 


4-6 06 19 

\^ kj kjkj J. <y 


13.0 


2891 3 

KJ . iUJ J-'-J 


36 


t;n(~ii"+ t^ilsti LiLi tii 1 1 1; atrir 

OllVJlL LJ L- 1 1U Ll 1} Lllo Ct L Wl 


V20 




18 


26 


42.3 


4-7 05 02 


11.2 


05681 


20 


short period pulscitor 


V21 




18 


26 


19.7 


+7 27 58 


12.2 


14068 

KJ . 1WUO 


10 


OllLJIL IJ L.l 1U LI 1} LI 1 o cl IU1 


V22 




18 


26 


54.2 


4-6 58 06 

1 u ijo kjkj 


11.8 


14862 


15 


Grinrt tiliti Lif i tii 1 1 t;n tnr 

OllLJIL IJ L.l 1U LI ULllDdlUl 


V23 


NCC 6633 54 

I'VJv ' ' KJ ~ J' J ' ' i 


18 


26 


44.9 


+6 24 20 


10.1 


0.07472 


4 


short period pulsS/tor 


V24 


NOC 6633 89 


18 


27 


33.4 


4-6 56 00 

\^ kj kjkj 


10.7 


06734 

KJ . KJKJ 1 U1 


10 


Qrinvt T1LST1 LlM TiiilQarnT - 
ollL/1 L IJ L.l IL/Li LJ UloCLLWl 


V25 




18 


27 


40.9 


4_7 08 34 


11.4 


05363 


18 


short period p ulscitor 


V26 


NCC 6633 108 


18 


28 


04.0 


4-5 58 13 


10.3 


08375 

KJ . UOO 1 'J 


5 


ehnrt Tif^TinH niiKafAV 

ollL/l L ULl Iw LI |J LI 1 o cL L W 1 


V27 


OSO 00445-00634 

v_TkJ\_/ KJKJ'-i'-t'J UUUO^ 


18 


27 


53.3 


4-6 08 51 

I LJ UO 'J A- 


11.2 


34196 

KJ . ij'-t L ,JKJ 


25 


KK T .VTUP 
LVLv LJ y 1 cTtL, 


V28 


CI* NCC 6633 VKP 188 

v^l 1 l VJT K^J KJKJ*-J*-J V IV! 1UU 


18 


28 


24.4 


+6 13 37 


11.3 


1.11068 


17 


eclipsing biriSyry 


V29 




18 


28 


57.3 


-l-fi 91 09 


13.3 


1.81161 


253 


LL^iijj&iiiti uiiioji y 


V30 


NOC 6633 1 31 


18 


28 


33.8 


4-6 53 1 5 

\^ KJ \JrJ J-lJ 


10.0 


04723 


8 


ollL/l L ULl Iw LI LJ LI 1 o cL L W 1 


V31 


NOC 6633 1 35 


18 


28 


42.9 


4-6 51 25 

|^ KJ \J -L iitj 


10.4 


03634 

KJ . UJUU^ 


(-; 


ennrt Tif^TinH nulQatni' 

ollL/l L ULl Iw LI LJ LI 1 o cl. L w 1 


V32 




18 


29 


10.2 


+6 43 52 


12.9 


0.11138 


24 


short period pulsator 


V33 




18 


28 


58.0 


+7 28 33 


12.0 


0.09568 


188 


short period pulsator 


V34 


NGC 6633 149 


18 


29 


17.3 


+5 39 20 


10.8 


0.04683 


6 


short period pulsator 


V35 


NGC 6633 147 (HD 170451) 


18 


29 


12.8 


+6 47 14 


9.5 


0.37532 


355 


W UMa 


V36 


V* BP Ser 


18 


30 


13.4 


+6 16 50 


12.5 






LB 


V37 




18 


30 


29.0 


+5 48 39 


12.8 


0.04203 


12 


short period pulsator 


V38 


GSC 00458-01442 


18 


30 


41.9 


+6 47 50 


10.9 


0.31777 


9 


RR Lyrae 


V39 


CSC 00458-00757 


18 


30 


59.8 


+7 11 51 


10.7 


0.11008 


1 


short period pulsator 


V40 




18 


31 


02.4 


+6 05 01 


12.3 


0.03941 


7 


short period pulsator 


V41 


V* BI Oph 


18 


31 


07.6 


+7 00 32 


11.2 


202.59 




Mira Cet 


V42 




18 


31 


18.7 


+7 05 25 


12.7 


0.07343 


9 


short period pulsator 


V43 




18 


31 


19.5 


+7 26 21 


13.3 


0.53548 


186 


RR Lyrae 


V44 




18 


31 


59.2 


+5 40 19 


10.7 


0.43263 


122 


RR Lyrae 


V45 




18 


32 


05.6 


+7 14 56 


11.9 


0.05936 


9 


short period pulsator 


V46 




18 


32 


51.9 


+6 49 01 


13.3 


0.50552 


156 


W UMa 


V47 




18 


33 


17.3 


+5 58 31 


12.2 


0.21377 


146 


pulsating 


V48 


GSC 00458-01108 


18 


33 


16.4 


+6 29 08 


9.1 


0.15686 


6 


short period pulsator 


V49 




18 


34 


03.5 


+6 36 03 


12.6 


0.13636 


12 


short period pulsator 



(v) neglecting the effect of blends in a crowded field; 



(vi) ignoring the possible effects of stellar companions on 
detection p robability (only ~ 1/3 of l ate- type dwarfs are 
single stars. IPuauennov fc Mavo'rlllQQ J) . 



Therefore this is only an order-of-magnitude figure, but will 
provide a useful comparison to the results of our more de- 
tailed analysis. 

We note that iHornel (l2003h estimates a detection rate 
of 3 planets per month for our project, using simple scaling 
laws. However, by his calculations, all the transit surveys 
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Figure 7. The probability of observing at least one (solid line) 
or two (dotted line) complete transits of a transiting planet in the 
lightcurves we have obtained for the NGC 6633 field. A minimum 
of 100 in-transit data points are required for a detection (so that 
10-mmag precision yields a total signal-to-noise ratio of ^ 10 for 
a typical transiting hot Jupiter). A typical transit duration of 
0.1 day is assumed. The one-transit probability also includes the 
requirement of at least one additional partial transit, so that the 
period can be estimated. The mean probability over this period 
range (1-9 days, sampled logarithmically) is 0.4 for one transit, 
and 0.3 for two transits. 

combined should be detecting 191 planets each month, which 
is clearly an overestimate. The main goal of these calcula- 
tions was perhaps not an estimate of the absolute detection 
rate, but rather a comparison of the various projects. If we 
assume the six planets found by the transit method to date 
were detected from a single year of observations by all the 
projects, and scale Home's estimate to match this rate, then 
his estimate for our search becomes 0.09 planets per year. 
For a transit search with para meters compar able to ours 
(though a shorter observing run) . iBrownl i2003l) estimates a 
rate of ~ 0.7 detected planet per 10 4 stars monitored (to a 
magnitude limit of R — 12). Our estimate above translates 
to ~ 2 planets per 10 4 stars. Brown's analysis includes a 
number of effects, such as that of the planet period and ra- 
dius distributions, which reduce the detection rate compared 
to our simple estimate. 



7.2 Monte Carlo simulation 

In order to obtain a more realistic estimate, we have mod- 
elled our observations using a Monte Carlo simulation. This 
allows us to use a realistic photometric noise model, and to 
take various observational effects into account (in particular 
blending, and the set of epochs at which the lightcurves are 
sampled) . 

We begin by simulating a population of stars in a 1 
square degree field of view. We use the luminosity fun ction 
for th e Solar neighbourhood, $o(A'/y-) (as tabulated in lCoxl 
l2000l . for My = —7, —6, ...18), allowing for an exponential 
drop-off with height above the G alactic plane (with scale 
heights, H(M V ), from lAllenlll975t . The number of stars of 
absolute magnitude Mv per cubic parsec, at a distance d 
and Galactic latitude b, is thus 



$(M v ,d,b) = $ (Mv)e" d|sinl|/fl(%) (6) 

We divide the volume sampled by a 1 square degree 
field into a series of shells of thickness 20 pc. We fill each 
shell with the appropriate number of stars in each My bin. 
To calculate the apparent (V) magnitude of each star, we 
assume an extinction of 1.9 magnitudes per kiloparsec in 
the Galactic plane, dropping off exponentially with distanc e 
from the plane, with a scale height of 140 pc iAllenlll973l) . 
Only stars with 5 < V < 19 are included in the model, as 
stars brighter than this can be avoided during field selection, 
while fainter stars are not resolved in APT images. 

Depending on Mv , an appropriate fraction of stars (tab- 
ulated in ICoxl boOO) are simulated as binary systems. For 
simplicity, a luminosity ratio is randomly selected from a 
uniform distribution in the range (0, l) 6 . The absolute mag- 
nitude of the primary is adjusted for this (keeping the total 
luminosity of the system equal to the original Mv), and the 
secondary is ignored. 

Apart from binary companions, a 3-pixel radius pho- 
tometry aperture centred on a star in an APT image typ- 
ically includes light from several other stars. We randomly 
assign pixel coordinates to the simulated stars, within a sec- 
tion of the APT CCD corresponding to one square degree. 
For each star, we then calculate what fraction of the total 
starlight in the aperture comes from that star, and take this 
into account when calculating the observed transit depth 
due to a planet. Stars for which this fraction is less than 5% 
are excluded from further analysis, as any planetary transit 
signal from these stars would be diluted below our detection 
limit. 

Around each remaining star, we place a planet with 
a radius of R p = 1.2 Rj up . We randomly select a period 
in the range 1-9 d ays, using the distribution described by 
iGaudi et al.l i2005l) : uniform in log P in each of two bins, 1- 
3 d ("very hot Jupiters", VHJ) and 3-9 d ("hot Jupiters", 
HJ), with the frequency of VHJ relative to HJ being 15%. 

Tables in iLand Jl992l) give the mass (M*) and radius 
(R») of the star based on its absolute magnitude, assuming 
it is on the main sequence. For the planet, we assume a circu- 
lar orbit, with radius a = (1 AU)(M»/M ) 1/3 (P/(1 yr)) 2/3 . 
We then randomly select an orbital inclination (i) leading 
to a transit (i.e. select cosi from a uniform distribution in 
the range 0-Pt r , where P tr ~ R*/a is the probability of such 
an inclination for a randomly oriented orbit). 

Neglecting limb-darkening, grazing transits, and 
ingress/egress times, we approximate the signal as a 
rectangular transit of depth AF = (R p /R*) 2 , and duration 

At = — a/(R„ +R p ) 2 - (a cosi) 2 (7) 
ira 

We assume that we will reliably detect planets with AF > 
AFmin = 10 mmag, provided multiple transits are ob- 
served and a minimum total signal-to-noise ratio is reached 
((SyW) m i n = 10). Considering that the lightcurves we have 



6 The observed distribution of binary mass ratios peaks 
near q = 0.3 and drops off considerably at higher values 
iDuQuennov fc M ayor 1991). Therefore on average, our simpli- 
fication overestimates the dilution of the primary's light by the 
companion, leading to a slightly lower estimate of the probability 
of detecting a planet around the primary. 
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obtained to date frequently contain systematic trends of or- 
der 10 mmag, this value of AF m in is somewhat optimistic. 
For our present analysis of the NGC 6633 lightcurves, a value 
of 20 mmag would be more realistic. Well-sampled transits 
of that depth (superimposed on the systematics) are eas- 
ily detectable, at least by visual inspection. However, we 
are co nfident that our implementation of the lKovacs et all 
J2004h trend filtering algorithm, combined with transit de- 
tection software, will allow us to reliably detect transits as 
shallow as 10 mmag (Sec. 15.41 . Thus we will continue to 
use this optimistic value for AF m i n , while noting that it is 
an important parameter. Its effect on the planet catch is 
discussed below. 

We estimate the photometric precision per measure- 
ment using a noise model which includes Poisson noise due 
to star and sky flux, and a limiting precision of 1.5 mmag 
for the bright stars (the solid curve in Fig. ■ Combined 
with AF, this tells us the number of in-transit measure- 
ments (iV m i n ) that need to be made in order to achieve a to- 
tal signal-to-noise ratio (S/N) m i u . Given a simulated planet 
and a set of observation epochs, we can calculate the prob- 
ability (P bs) that a minimum of 2 complete transits will be 
observed, including at least AT m i n in-transit points. We do 
this by testing 1000 transit epochs. 

To estimate the total planet catch for the assumed 
observing parameters (Galactic latitude and observation 
epochs), we add up the values of Ptr x Pobs for each sim- 
ulated planet, and scale the result to the appropriate field 
size. We also multiply by the fraction o f stars that do hos t 
short-period planets, which is ~ 1% (e.g. lGaudi et all2005l) . 
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Figure 8. Distribution of period (top panel), host star abso- 
lute magnitude (middle panel) and apparent magnitude (bottom 
panel) for detected extrasolar planets in the Monte Carlo simula- 
tion. The histograms were generated by weighting each simulated 
planet by the probability of detecting it (P tr x Pobs)- The sim- 
ulation parameters were set to match those of our observations 
of the NGC 6633 field in 2002-2004. In each panel wc show the 
distribution for three different values of the minimum detectable 
transit depth. 



7.3 Validating the simulation 

In order to ensure that this rather complicated simulation 
yields reasonable results, we performed a number of tests. 
Firstly, we generated a set of fake APT images using the 
distribution of stars from the simulation (for b « 8°). The 
simulated images appear very similar to a real image at the 
same Galactic latitude (with the exception of some patches 
of increased extinction in the real image). Processing these 
fake images through the photometry pipeline, we also con- 
firmed that the star counts as a function of limiting magni- 
tude closely match those obtained from a real image. 

We then modified the si mulation to rep roduce the sim- 
plifying assumptions made bv lHornel i2003h . The result was 
within a factor of 2 of Home's estimate for our project. Fi- 
nally, we confirmed that the result scales with planet, star, 
and observational parameters (R p , a, My, R*, sky bright- 
ness, survey duration and (S/N)min) as predicted by Home's 
equation 3. 



7.4 Results for the NGC 6633 field 

We have simulated observations at the Galactic latitude of 
our NGC 6633 field (b — 8°), using the epochs of the best 
images we have obtained for this field in 2002-2004 (2486 
images, from 61 nights over 3 seasons). The expected planet 
catch for a single 6-square-degree APT field is ~ 0.1 planet 
(AF m i n = 10 mmag, 2 complete transits). If we accept one 
complete observed transit as a detection (plus at least one 
partial transit so that some estimate of the period can be 



made), the estimate becomes ~ 0.2 planet, which agrees well 
with the rough estimate derived in Section f7. II 

As expected, we are most sensitive to planets with the 
shortest periods (Fig.|HJ. However, since we have simulated 
planets with periods longer than 3 d to be almost 10 times 
more frequent than those with P < 3 d (according to the 
observed distribution), the distribution of detected planets 
also has a peak at 3-5 d. Most detected planet hosts have an 
apparent magnitude V w 13, while almost none are fainter 
than V w 14.5. 

The luminosity distribution of detected host stars peaks 
at My ~ 4 (for AFmin = 10 m mag), as also predicted by 
iPepper. Gould fc Depovl i2005t) . Amongst the stars that are 
small enough for a 1.2 Rj up planet to cause a 10 mmag 
transit, those with My « 4 are the brightest and there- 
fore the most numerous in our magnitude-limited sample. 
As a result, AF m - m is perhaps the most important parame- 
ter in determining the total planet catch. 7 In the case of the 
NGC 6633 dataset, about 90% of the systems detected with 
AFmin = 10 mmag have transits shallower than 20 mmag 
(Table |2J . This highlights the importance of removing sys- 
tematic trends from our lightcurves so that shallow signals 
can be detected. 

Repeating the simulation for one of the other 7 fields 
(for each of which we have ~ 1000 useful images, over 2 

7 The largest detectable host-star radius is AF m i n — 1//2 R P . Vary- 
ing the assumed planet size (Rp) is equivalent to varying 
AF ■ - 1 / 2 
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Table 5. The total number of detected planets per APT field 
(minimum two complete observed transits) as a function of the 
minimum detectable transit depth (AF m ; n , given in the first row 
in mmag). We compare two simulations: our NGC 6633 observa- 
tions (2002-2004, V band, Sec. 16.11 . and higher-cadence obser- 
vations at a Galactic latitude of 25° in / band over two months 
(Sec. 1731. 
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seasons) gives 0.03 planets per field (AKnin = 10 mmag, 2 
complete transits). The reduced estimate compared to the 
NGC 6633 field is due to the smaller number of images, 
fewer nights with long continuous coverage (mostly due to 
bad weather), and sparser time sampling. Thus our most op- 
timistic estimate of the planet catch in the data we collected 
over 2002-2004 is ~ 0.3 planets, or ~ 0.6 if we settle for one 
complete transit. 




2 0.06 - 
-o 

% 04 ~ 
a) 

S 0.02 - 




0.6 

log(P) (days) 



- 1 


n 




1 


i i i i 1 i i i - 
&F m{n - 10 mmag 










AF = 15 mmag — 










AF - ?n mmao " 








1 


' — l_ 

— i — i — 



7.5 Discussion 

The above disappointing estimate suggests that our initial 
strategy of simultaneously monitoring four fields over mul- 
tiple seasons is far from optimal. We selected this strategy 
believing it would maximise the number of stars we monitor, 
while obtaining a minimum time-sampling required to detect 
a transit (one image every ~ 14 minutes). The signal-to- 
noise ratio (S/N) of a detection would be increased by fold- 
ing the lightcurve at the appropriate period. However, gaps 
in the observations due to weather and technical problems 
have reduced the probability of detecting multiple, complete 
transits (Fig. UJ . To date, the systematic trends present in 
many lightcurves have prevented the use of detection soft- 
ware to make shallow transits detectable by phase-folding. 
We need to sample the lightcurves more frequently so that 
the required total S/N can be accumulated over a small 
number of observed transits. 

We have used the Monte Carlo simulations to study 
the effects of various observational parameters on the planet 
catch, and thus select a more efficient search strategy. One 
important change was the choice of filter. Other transit- 
search teams have discussed the merits of observing at long 
wave lengths, in particular in I band (e.g. Ivon Braun et all 
2005). In terms of maximising the detection rate, there are 
two main advantages. First, the greater CCD sensitivity in I 
band allows shorter exposures (60 s, compared to 150 s in V) 
and thus higher-cadence observations. Second, observing in 
I band significantly increases the number of well-measured 
red dwarf stars in our sample (which have deeper transits) 
relative to their larger, bluer counterparts. This makes the 
total detection rate less sensitive to the minimum detectable 
transit depth (see Tableland Fig. [§J|. 

The simulations indicate that the detection rate is es- 
sentially constant from the Galactic plane out to a latitude 
of |6| w 25°, and is only reduced by ~ 20% even at |6| « 65°. 
This is due to the reduced blending in less crowded fields, as 
well as the near-isotropic distribution of the well-measured 
red stars (most of which are within ~ 100 pc). Although we 
have not calculated the frequency of "false positive" detec- 



Figure 9. Distribution of the properties of detected systems for 
a model of our optimal observing strategy. Compare to Fig. [5] 
representing our initial strategy. We are confident that with future 
improvements in our trend filtering and detection software, the 
AF m i n = 10 mmag case will become a realistic one. 



tions due to eclipsing binaries blended with a third star along 
the line of sight, this will be significantly lower in the higher 
latitude fields. Furthermore, the photometric precision (and 
absence of systematics) we have assumed in the simulation 
are easier to achieve in less crowded fields. For these reasons 
we conclude that fields in the range 15° < \b\ < 45° are 
optimal. 

Future fields will be at southern declinations, so that 
they are observable at low airmass for a full 8 hours every 
night. 8 We will select pairs of adjacent fields, obtaining an 
image of each field every 4 minutes. 

We have simulated this new strategy, taking into ac- 
count the different stellar magnitudes, lower interstellar ex- 
tinction and brighter sky in I band. We can achieve a de- 
tection rate of 0.2 planet per field (A_F m i n = 10 mmag, 2 
complete transits) in only 20 nights of observations (taken 
over up to 2 months, allowing for gaps due to bad weather 
and full Moon). 

Observing 10 fields per year (one pair for every 2 
months, allowing for telescope downtime), will lead to an 
overall detection rate of 2 planets per year, a significant 
improvement over the estimated 0.2-0.3 per year for our 
2002-2004 observations. If further efforts do not allow us 
to reliably detect transits shallower than 20 mmag, the es- 
timate for the new strategy is halved, while for the past 
observations it is reduced by a factor of 10 (Table Q>J. 



8 From SSO, we have a minimum of 8 hours of darkness per night 
for 8 months each year. 
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8 CONCLUSIONS 

We have described in detail the methods we have used to 
achieve millimagnitude-precision relative photometry using 
the Automated Patrol Telescope. For bright stars, the pre- 
cision was limited by intra-pixel variations in CCD sensi- 
tivity (IPSV), coupled with undersampling of the images. 
This was largely overcome by a new observing technique, 
which optimally alters the effective PSF in order to average 
out the effects of IPSV. We have demonstrated that using 
this technique, a customised aperture photometry package 
and an automated processing pipeline, we can obtain the 
large number of high-quality lightcurves required to detect 
transiting extrasolar planets. 

During 2002-2004, we observed 8 fields near the Galac- 
tic plane. Our analysis of data from the field centred on the 
open cluster NGC 6633 yielded 4 transit candidates. Follow- 
up observations revealed each of these to be eclipsing binary 
systems. We also found 49 variable stars in this field. Three 
of the other fields also contained binary stars with eclipses 
similar to those of a transiting planet. 

By simulating the population of stars we observe, the 
properties of the telescope, and our observing strategy, we 
have studied the sensitivity of our search to short-period 
planets. Assuming the frequency of planets around stars 
targeted by radial velocity surveys is representative of the 
Galactic field stars we observe, we have estimated our 
expected detection rate. Considering all the good-quality 
lightcurves for the 8 fields we have observed to date, we 
expect to find 0.3-0.6 planets. 

Our new strategy, whereby we obtain high-cadence ob- 
servations in I band, will yield an estimated 2 planets per 
year, a 1 0-fold increase in detection rate over previous obser- 
vations. We have been using this strategy since November 
2004. The new camera for the APT, being built in 2005, 
will have at least 3 times the field of view of the current 
one, higher sensitivity, and smaller pixels. This will further 
increase the detection rate by a factor of a few. 
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